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Abstract 

We consider a system made up of electrons interacting with a neutralizing positive background 
within a cubic box of volume V . After dividing the box into N (or cubic cells for the polarized 
(unpolarized) case, we average the creation field operator over each cell with a suitable weight 
function and we consider the quantum crystalline states obtained by letting all the average operators 
act on the vacuum state. These states exclude the possibility that each cell may momentarily 
contain more than one or two electrons in the polarized or unpolarized case. The expectation value 
of the Hamiltonian over this class of states is evaluated in the thermodynamic limit and the weight 
function is chosen in such a way to minimize the expectation value. The involved numerical analysis 
is explicitly performed with a weight function having a generalized Gaussian shape depending on 
a parameter. It turns out that the unpolarized and polarized quantum crystalline states yield an 
energy per particle smaller than the homogeneous Hartree-Fock ones for > 90 and Vg > 28, 
respectively. Moreover, for the polarized case, the energy per particle at = 100 is -0.01448ryd 
close to -0.0153530(8)ryd, the best quantum Monte Carlo value [Drummond et al, Phys. Rev.B 
69, 085116, (2004)] and this discrepancy measures the correlation contribution neglected in our 
approximation. 
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I. INTRODUCTION 



Quantum Coulombian systems determine the properties of common matter but their 
theoretical investigation is far from being complete. Around the sixties it has been proved 
that they are H-stable^'^ and thermodynamically stable^ provided that they are overall 
neutral and all the constituting species with electrical charges of a given sign are made up of 
Fermions^'^. From this result follows that the so-called jellium model^ of metallic conductors 
has a fundamental state |\E'o) with eigenvalue Eo{N, V) such that eo(n) = Eo{N, V)/N exists 
as function of the particle number density n = N/V in the thermodynamic limit: N ^ oo, 
y — > oo with n fixed. Choosing units such that h = 1, the Hamiltonian of the model reads 

h = 4 + ^ = EE|^<^«m+ (1) 

k s=l 

1 X - 47re^ I I 

k,q,p ^ ^«k+q,.«p-q,s'«P,«'«k,. 
s,s'=l ^ 

where the involved symbols have the standard meaning (see, e.g., section 3 of Ref. p]) and 
the prime on the summation symbol means that the value q = is excluded from the sum, 
a convention adopted throughout the paper. Wigner^ evaluated the expectation value of H 
over the state 

iF„)= n ri<.io)' (2) 

|k|<fcp ^=1 

describing the fundamental state of a degenerate Fermi gas of (non-interacting) electrons 
with Fermi momentum = (Svr^n)^/^. The resulting energy expression coincides with that 
obtained the Hartree-Fock (HF) equations under the homegeneity assumption. It reads 

e,,„(r,) = (F„|H|F„)/(iVeV2aJ (3) 

2.21 0.916 
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where 

rs = rjag, ao = l/me^, 
r„ = (3/47^)1/3 _ (97r/4)i/7A:p, (4) 

Ts being the perturbative parameter and m the electron mass. On general grounds one has 

^o{rs) = €-Jl,u{rs) + €-coTr,u{rs), (5) 
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where ecorr,u{fs) denotes the sum of the remaining terms in the perturbative expansion. It 
obeys the inequahty ecorr,u{rs) < and is named correlation energy. Its evaluation is by no 
way easy. In fact, no further progress in its analytical knowledge was made since papers ''"^^ 
that yielded the expression 



the last addend being the most recently evaluated one^^. Expression (6) can only be accu- 
rate for dense systems, because at low density it becomes positive and violates the reported 
inequality. Even though €cott{i^s) is practically unknown from an analytical point of view, 
it is expected to be small in comparison to Eq. (3). In fact, on the one hand the min- 
imum of Eq. (3) and the corresponding value compare favourably with the ionization 
energy and the value of metallic sodium^. The same happens for the bulk moduli and 
the cohesive energies of some typical metals (see table 2.2 of Ref. [^^] and table 5.3 of [^^]). 
On the other hand, density functional theory (DFT)^^ and quantum Monte Carlo (QMC) 
calculations^^"^^ nowadays provide accurate numerical approximations of the correlation en- 
ergy that definitely results to be small in comparison to Eq. (3) [see Figs. 1 and 2]. This 
property also applies to the polarized jellium. In fact, Bloch^^ showed that, if the funda- 
mental state of the jellium is assumed to be fully polarized, the simple HF approximation 
of the energy per particle is 



This energy becomes smaller than Eq. (3)'s as Vg > 5.7. Thus one expects that the jellium, 
as its density decreases, passes from the unpolarized to the polarized state. Moreover, at 
low densities, both Eq. (3) and Eq. (7) yield an energy per particle greater than that of ionic 
crystals. This observation led Wigner^ to suggest that, at low density, the jellium becomes 
a crystal (known as Wigner crystal) in the sense that each of its electrons oscillates around 
its equilibrium position, forming in this way a lattice of harmonic oscillators. In fact, under 
the assumptions that the crystalline structure is a bcc one and that the harmonic oscillators 
are decoupled, Wigner^^ obtained the following energy per particle 



ecorr,u{^s) = 0.0622 lu r, - 0.094 + 0.018 r, Inr, 
- 0.020 r, + 0(r,2 In r,). 



(6) 




2^/3 2.21 21/^0.916 
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Ts > 1. 
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FIG. 1: The energy per particle of the jelhum model (in Rydberg units) in the inner part of the 
r^-range according to different approximations. Unpolarized case: Eq. (3) — > dotted curve, QMC 
results^^ full circles; SLQCS long-dashed curve. Polarized case: Eq. (7) — > dot-short-dash; 
Eq. (9) —> dot-long-dash; QMC results^^ ^ full triangles: SLQCS— > continuous curve. 

Carr^^ improved this analysis taking into account the coupling among the oscillators and 
(leaving aside a further small positive 0(r7^) term) found that 

^ 2 / N 1-79 2.66 

^crst/(iVeV2aJ = + rs > 1. (9) 

rs rs 

Both expressions are definitely smaller than Eq.s (3) and (7) at large (see Figs. 1 and 
2). For a reason explained below, this finding would not have been sufficient to conclude that 
Wigner's suggestion holds true. This conclusion only follows from the fact that DFT^^ and 
QMC analyses^'''^^'^''"^^ definitely showed that, at high dilution (i.e. r^, > 106±1 according to 
the most recent QMC analysis^^), the stable phase of the jellium model is the fully polarized 
bcc one. Hence, Wigner's idea that the low density jellium model has a crystalline structure 
presently can be considered fully established. In passing we mention that this idea has 
recently found further applications in relating two-dimensional Wigner crystals in presence 
of a strong magnetic field to quantum Hall effects^^"^^. The reason why Wigner's conclusion 
cannot be drawn from the fact that the energies predicted by Eqs. (8) and (9) are smaller 
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FIG. 2: Outer r^-range. The curves have the same meaning as in Fig.l. The full squares are the 
QMC values^^ relevant to the polarized bbc phase. (One full square and all full triangles are hardly 
distinguishable from the corresponding full circles.) 

than that of Eq. (7) is related to the fact that, in deriving Eqs. (8) or (9), Wigner and Carr 
respectively substituted Hamiltonian (1) with that of appropriate uncoupled or coupled 
oscillators plus an appropriate constant. By so doing, the Hamiltonians do not coincide with 
Eq. (1) and one no longer can invoke Rayleigh-Ritz' principle to conclude that the quantum 
state associated to the crystalline picture is closer than Eq. (7)'s to the true fundamental 
state. This appears even more evident from the fact that the energies predicted by Eqs. (8) 
and (9) are smaller than Eq. (7)'s throughout the full range. In principle, to show that 
at low density the fundamental state of (1) corresponds to that of a crystal, one should 
show the the minimum expectation value of H only occurs over an appropriate normalized 
quantum crystalline state |\l/cr,o)- Of course, this would correspond to exactly solve the 
problem. This being - if not impossible - very hard, one has to confine himself to look for 
a quantum crystalline state such that the expectation value of H on it is smaller than that 
over state |Fu). Actually, this has already been done in Ref. p'^] by solving the HF equations 
for a Wigner crystal due to the fact that the HF equations coincide with equations obtained 
minimizing the expectation value of H over a determinantal wave-function^^. In this paper 
we shall provide a different procedure that, based on Rayleigh-Ritz's principle and quantum 
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field theory, allows us to determine a quantum crystalline state with an energy smaller than 
homogeneous HF's. To this aim, in section IIA, we shall first average the creation field 
operator ipa{r) over each cell of the crystal with a weight function i/(r) having the crystal 
periodicity. Denoting the resulting operators hj tpj^^, these obey the Jordan- Wigner algebra 
[see Eq.s (13), (16) and (17)]. Hence, the Rayleigh-Ritz principle is appplied to the class 
of normalized antisymmetric quantum cristalline states |\l/cr) = Ilma'^omlO) determine 
the unknown z/(r) that, in turn, determines the electron density function n(r). These cr)s 
will be named strongly localized quantum crystalline states (SLQCS) because they neglect 
the possibility that two or more electrons may momentarily be present in the same cell. 
Then, the expectation value of the kinetic and potential operators are evaluated over the 
SLQCS in terms of the electron density in sects. II. B and II. C for the unpolarized case and 
in sect. Ill for the polarized one. The requirement that the resulting expectation value of H 
be minimum determine the equations that determine n(r). In sect. IV we show that these 
equations coincide with the HF ones with a different boundary condition and this implies 
that exchange and correlation contributions are equal to zero so that electron interactions 
are not fully accounted for. Aiming to avoid, as far as possible, computer calculations we 
do not attempt to solve the mentioned equations. We simply consider a family of electron 
density functions, depending on a parameter a, and numerically evaluate the expectation 
value of H by the formulae derived in sects. II and III over a grid of a values. Then, for 
each Ts value, we look for the a value which makes the expectation value of H minimum. 
The results of this analysis are partly shown in Figs. 1 and 2 and discussed in more detail 
in sect. V through their comparison with HF and QMC values. Our final conclusions are 
drawn in sect. VI. 

II. BASIC RELATIONS 

A. Strongly localized quantum crystalline states 

In order to construct a quantum state that describes the electrons in a crystalline config- 
uration we shall proceed as follows. At first we shall assume that the electrons are confined 
within a cubic box (centered at the origin of a Cartesian frame) of edge L and volume 
V = L^, and that the system is not polarized. The cube is divided into cubic cells of 
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edge a where Mc denotes the number of cells along Cartesian axes x, y and z (or Xi, with 
z = 1, 2, 3). To relate the values of Mc and a to and V we need some preliminary remarks. 
Each cell generally contains two electrons, respectively with spin up and down, in order to 
fulfill the unpolarization condition. This condition can exactly be fulfilled only if N/2 = M^. 
It is however clear that this restriction fades away in the thermodynamic limit. In fact, one 
can generally write = 2Mj} + A'^,. where Mc is the largest positive integer such that the 
previous relation holds true with Nr > and integer. According to the above relation, 2M^ 
electrons are accommodated in the M^ cells that together form the outset cubic box V, 
while the remaining A^,. electrons will be located in Nr/2 or Nr/2 + 1 (depending on whether 
Nr is even or odd) cells next to the external boundary of V. Recalling that the electric total 
charge of each cell (but one, if A^^ is odd) is equal to zero, the contribution of these border 
cells becomes negligible in the thermodynamic limit (see pages 402-404 of Ref. [^^])- As 
A^ ^ oo we have that Mc — ^ oo and from the above considerations follows that the lattice 
spacing is 

a = L/Mc ^ L{2/Nf''^ = (87r/3)^/Vo. (10) 

For notational simplicity, in the following we explicitly refers to the case of odd Mc so as to 
write Mc = 2M + 1. Then, Vm, the m[= (mi, m2, m3)]th cell of V, is defined as 

Vm ={r I (mi - ^)a < Xi < {nii + ^)a, 

mi = -M,...,M, i = 1,2,3}. (11) 

By construction, 

V= U V^= [j (12) 

-M< 

where Z^y = {~\V is the intersection of the cubic three-dimensional lattice Z^ of spacing 
a with the outset box V . Introduce now the operators 

V'a,m = I i^ir)^air)dv, ipl^^ = / V{r)4jl{v)dv, (13) 

obtained averaging the electron field operators of destruction and creation over the mth cell 
respectively with a weight function z/(r) and its complex conjugate ^'(r). Function v{v) is 
assumed to be periodic with period a, i.e. 

z/(r) = i^(r + am), (14) 



7 



m being an arbitrary triple of relative integers. Moreover we also require that i/(r) and all 
its first partial derivatives are continuous throughout the closed unit cell and, finally, that 
i^(r) has unit norm over the unit cell, i.e. 

I \v{r)\'^dv = 1. (15) 

Definitions (13) and the canonical anticommutation relation of fields ifja and ipl^ yield the 
following anticommutation relations 

{^«,m, V'^^m'} = ^a,/3'^m,m', (16) 

One immediately verifies that the state "V^^ ^nlO) contains one electron within the mth cell with 
spin up or down depending on whether a equals one or two. Furthermore, the expectation 
value of the particle number density operator n(r) over the previous state is 

{0\^Pa,m nir) <^|0) = |z/(r)pe(reV;,) 

= n{r)Q{reV^), (18) 

where, similarly to Heaviside's function definition , 6(r G Kn) is defined to be equal to one 
or zero depending on whether the tip of r lies within or outside Vm- Eq- (18) shows that 
?7,(r), the square modulus of function z/(r) entering definition (13), times 0(r G Kn) is the 
particle number density of the electrons within the full box for the quantum state "ipa^lO). 
Consider now the fully antisymmetric state 

iwc.)=( n ti^u)\o) (19) 

—M<mi,m2,m[i<M a=l 

and label the index pair (a, m) by a single index J' that runs over {1, . . . , A^}. Then Eq. (19) 
is more compactly written as 

N 

|Wcr)= (n^:7)|0)- (20) 

J=l 

One easily verifies that (WCT|?^(r) |Wcr) = n{r) [the latter's value on the cells' boundaries 
being defined by continuity]. Hence, quantum state |Wcr) describes a unpolarized crystalline 
configuration of the electrons within the cubic box V with a periodic density n(r) = 
|i^(r)p. It should be noted that the states of form (19) have a peculiar feature: they do 
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not allow that two (or more) creation field operators, with the same spin, act on the same 
cell. This limitation will be discussed in more detail at the end of section IV. As already 
anticipated, quantum states of form (20) are named strongly localized quantum crystalline 
states (SLQCS). Despite this fact, it looks reasonable that at high dilution the fundamental 
state of the jellium model is closer to a SLQCS than to Eq. (2). Our task now is to find 
function z/(r) and the parameter region where the above property is fulfilled. This will 
be accomplished by: i) evaluating the expectation value of (1) over |Wcr), ii) performing 
the thermodynamic limit, iii) deriving the equations that determine the i^(r) function that 
minimizes the expectation value and iv) approximately solving the last equations so as to 
make the comparison of the resulting energy with homogeneous HF's possible. 

The evaluation of (Wcr|H|Wcr) requires that of (Wci-|T|Wcr) and (Wcr|V^|Wcr), where 



B. Evaluation of the expectation value of T 

The evaluation of (Wcr|T|Wcr) is easily performed by the following anticommutation 
relations 




(21) 



and 




k,p,q s,s'=l 



(22) 



a, (3 = 1,2, meZ^ 



(23) 



and 



{^^(r),^^,^} = 5,,^i^(r)e(r G V^), 
a, 13 =1,2, mGZ^ 



(24) 



the remaining anticommutators being equal to zero, i.e. 



{^a(r),V^W = {^i(r),^t_^} = 0. 



(25) 
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By the relabel m) J the above anticommutators are more compactly written as 

{V'a(r),^/';^} =F,,^(r) = b^^(iU{j)Q{j G V^), (26) 

with JT" = 1, . . . , and a = 1, 2. The subsequent use of Eqs. (26) and (20) yields 

N 

(27) 

Its adjoint determines (Wcr|'?Aa(r). In this way, by Eqs. (16), one finds that 

J,J'=l 



2M 
2M 

{o\ipN ■ ■ ■ ■ ■ ■ ■ ■ ■ V'^r-i^^y+i ■ ■ ■ ^]vlo) 

^ VF.,^(r) ■ VF..^(r) 

{oitpN ■ ■ ■ ■ ■ ■ ■ ■ ■ ^aJvIo) 

_ Vi^a.^(r) ■ Vi^.,^(r) _ ^ 
=1 2M 

(28) 

In obtaining the last relation, we converted JT" to (/9, m) and used Eq. (26). Finally, in 
the thermodynamic limit, the periodicity of z/(r) yields the sought for expression of the T 
expectation value 

TV r 

(W„|T|W„) = ^J^ \Vuir)\'dv, (29) 

where Vq is the cell defined by Eq. (11) with m = 0. It is now convenient to introduce the 
following dimensionless quantities 

f=2r/a, and u,{0 = ici/2f^^u{a^/2) = {a/2f/Mr). (30) 

^ is a dimensionless position vector such that the largest and smallest value of each of its 
components respectively are +1 and -1 as r ranges over cell Vq. In this way, ^ is confined to 
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vary within the cubic cell Vq (different from Vq) defined as 

Vo = {fl - l<e. < l,z = 1,2,3}. (31) 

By Eqs. (30) and (15) it results that i^i(0 has unit norm over Vq. Further, i^i(^) can be 
defined throughout imposing the periodicity condition i^ii^) = + 2m). In terms of 
the just defined dimensionless quantities Eq. (29) becomes 

(W„|T|W,r) = N^^(-f' [ \VMO\'d'^ (32) 

by Eqs. (32) and (10). 



C. Evaluation of the expectation value of V 

We find it advantageous to start from expression (22) of V because this accounts for the 
contribution of the positive background through the fact that the mode q = is excluded 
in the first sum of Eq. (22). Using the decomposition 

^a(r) = -i^(iKsVa, (33) 

k s=l 

from Eq. (13) one obtains 

^".m = EE^ / ^^'■Mr)rf^;. (34) 

The Fourier transform (FT) of restricted to the unit cell centered at the origin, is 

defined as 

i>(k) ~ [ iy{r)e-'^-''dv, (35) 

Jvo 

and obeys 



j){k) = z/(-k) (36) 

where, hereafter, the tilde denotes the FT and the overbar the complex conjugate. The 
periodicity of z/(r) allows us to write the integral present in Eq. (34) as e*°''' '"z/(— k), so that 
ipa,m reads 

'>Pa,m = V'J = ^J-XsO'Ks^ (37) 
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where we have put 

^J;Ks = • (38) 



From Eq. (37), its adjoint and the canonical anticommutation relations one finds that 

{a\,s,^j} = ^j;k,s and {ak,s,^j} = ^j-k,s- (39) 
By these relations the evaluation of (Wcr|V^|Wcr) becomes straightforward. In fact, 

N 

OvIWcr) = ^(-l)^-^^/'I . . . ^]j-l^J;Ks^ 

J=l 

^ki---^!v|0) (40) 

and 

N J-1 
J=l C=l 

ipl-.. ^LiV'i+i • • • '4'j-i'^j+i ■ ■ ■ ^ivlo) 

N N 
J=l C=J+1 

^i...^:^_i^:^+i...^Li^i+i...^]vio), (41) 

where it is understood that appended index values can neither be smaller than one nor 
exceed N . The above expression can be recast in the form 

N J-l 

'ak,.|We,) = 5^$^(-l)^+^(^^;k,.^z:;p,.' 

J=l C=l 



M=l 

and the prime over the product means that index M. never takes values J and C By 
relation (42) and its adjoint one gets 

N J'-l 

(Wcr|a^k+q,sa^p-q,s'ap,s'"k,s|Wcr) = X 

j>,j=i a=i 

i;(-l)^'^^'^^^^(v'^';k+q,.'V^z:';p-q,.' (43) 

C=l 

- '/^£';k+q,sV5j';p-q,s') (^J;k,s^£;p,s' " ^£;k,s^J;p,s') X 
(Wcr I^AT ■ ■ ■ ^J'+l^J'-l ■ ■ ■ 
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The matrix elements are different from zero only ii J' = J and C = C and, in this case, 
they are equal to one. Thus, one finds that 

( Wcr I a^k+q,sa'^p-q,s'ap,s"^k,s I Wcr) 
N J-l 

J=l C=l 

(v^J;k+q,s'V^£;p-q,s' ~ V'£;k+q,s</'j';p-q,s') ■ (44) 

The terms in the above sums are symmetric with respect to the exchange J ^ C and equal 
to zero if ^7 = £. Thus, we can let C range over 1, . . . , provided that we divide the result 
by two. Converting J and C to index pairs (a,m) and (a',m'), using definition (38) and 
performing the sums over a and a' one finds that 



_ u{-\^ - q);/(q - p);>(-k) ;/(-p) 

- Ks'Ks' (e-('^-P+^)-(— ') + C.C.)] , (45) 

where c.c. stands for complex conjugate. Each sum present in the above expression is the 
product of three sums of the form 

M M 



m=-M m'=-M 



One easily verifies that 



M M 



J2 J2 e^"''^"') = (2M + l)x 



m=-M m'=-M 

2M 2M 



p=-2M p=-2M 



In the limit M — > oo one finds that 

MM oo 

J2 e''"^^"'^ ^ (2M + 1) ^ e 



laqp 

m=—Mm'=—M p=—oo 



(2M + l)27r Y S{aq-2rmT). (48) 
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In the three dimensional case, recalhng that = = (2M + 1), one has 

J2 e»aq.(m^mO ^ ^{2^^ ^ ^(Qq_27rm). (49) 

Consequently, Eq. (45) becomes 

(Wcr |a"^k+q,sa"^p-q,s'ap,s'ak,s I Wcr) 

{2TifNi){-\i - q)i>(q - p)z>(-k) z>(-p) 

^ [5s,s5s,s5(aq - 27rm) (50) 

- 5s,s'5s,s'5(a(k + q - p) - 27rm)] . 
By this result, in the limit V oo the expectation value of V becomes 

(Wcr|V|W„) = ^x (51) 

[4 ^ '^(aq - 27rm) - 2 ^ 5(a(k + q - p) - 27rm) 

where value m = is excluded in the first sum because value q = is not allowed in 
Eq. (22), and the numerical factors in front of the sums account for the sum over s and s'. 
After introducing the dimensionless momenta 

P = ap/2, K = ak/2, and Q = aq/2, (52) 

from Eqs. (35) and (30) one obtains that 

MK)= [ e-''^<u,iC)d'^=(-fi>ik). (53) 

JVo ^'^^ 

After putting 

ni(a = ki(e)P = (a/2)Mr), (54) 

where the last equality follows from Eq. (30) and the fact that ?T,(r) = |i^(r)p, one easily 
shows that the following relations 

/ ^>i(K + P)z>i(K)t/3j^=(27r)3 f e*^-«>i(Orrf=^e 

J Jvo 

= (27r)3rri(-P) (55) 
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and 



u{P + 7rm)i>(P - Q + 7rm) (56) 



m 



£Z3 



8 / Wi{0\'e-''^-^^d'^ = 8ni{Q) 

JVn 



'Vo 

hold true. By Eqs. (52), (55), (56) and (10), Eq. (51) converts into 



(We.|\/|W,r) = ^-f-) X (57) 



2ao^s 7r Vtt 



2 , „ I ~ |2 



r ' pi(vrm) 1 f \ni{Q)\ 



Define the auto-correlation function of ni{^) as 



IQI 



9i{0= I ni{i+f])nMdS- (58) 

This differs from zero only within the cubic cell centered at the origin and with edge length 
equal to four. This cell will be denoted by 2Vq. One finds that 



'2Vb 

|2 



f e-^«-%(Orf='e=Ki(Q)r- (59) 

J2Vn 



Then, recalling that 47r/|Q| is the FT of 1/r, the integral in Eq. (57) can be written as 

|[-l(Q)]%3^_0.2 f 9l{0, 

Finally, Eq. (57) reads 



IO|2 d^Q = ^^' (60) 

1^1 J2}^ \Q 



Ne^ 1 / 3 \ 1/3 

(W„|V^|W„) = ;J^-(-) X (61) 



mG2 



2aors 7r Vtt/ 



III. FINAL EXPRESSIONS FOR THE UNPOLARIZED AND POLARIZED CASE 

In this section we report the final expressions of the expectation value of the Hamiltonian 
over the quantum crystalline states of the form (20) for the unpolarized and polarized cases. 
The first case has been analyzed in the previous section. There we found that, if we put 

th] ^ / \^^v^{i)fd% (62) 
15 



"*i-E'iiiiP (63) 



and 



VeH - - / '-^d% (64) 

J 2Vo 141 



the expectation value of H over the unpolarized crystaUine state defined by Eq. (20) is 

=2 



2 

e„[z/i,r,] = (Wer|H|W,r)„/ (iV^) (65) 



r 2 



with 



3\| l/3\^ /3 



In fact, the first contribution on the right hand side (rhs) of Eq. (65) represents the expec- 
tation value of T [see Eq. (32)], while the sum of the remaining two terms is that of V [see 
Eq. (61)]. It is stressed that expressions (62), (63) and (64) do not explicitly depend on the 
density of the system because they only involve dimensionless quantities i^i(0 gi{Cj. 
An implicit dependence is however present in ^'^(^i consequently, in ni{^), gi{C,) and 

^(K) so as to make the dependence of (\l/cr|H|\l/cr) different from that of Eq. (3). 

If we assume that the system is fully polarized, each cell exactly contains one electron 
with, say, spin up. Thus, the number of the cells is determined by the condition N = MJ^+N^ 
and the lattice constant is 

ap = L/M, ^ L/N'/' = (47r/3)i/Vo = (67) 

where a„ denotes now the unpolarized spacing, denoted by a in Eq. (10). The fully polarized 
and completely antisymmetric SLQCS becomes 

|Wcr,p)= ( n ^U)|0) (68) 

with Z^y denoting now the intersection of V with the cubic lattice of spacing a^. The 
calculation of (Wcr,p|T|Wcr,p) and (Wcr,p|V^|Wcr,p) proceeds similarly to what we did in 
sects. II. B and II. C with the difference that relabel (a,m) J becomes now (l,m) J . 
In Eq. (28), this implies that the sum over /5 no longer is present since we always have /3 = 1 
and that the sum over m yields instead of A^/2. Hence, the numerical factor on the rhs of 
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Eq. (29) is unchanged. However, when this integral is converted to the relevant dimensionless 
quantities, the factor a^t"^ converts into ap~^. Thus, recalling Eq. (67), numerical factor ct,u 
in the rhs of Eq. (66) must be multiplied by 2^/^. In evaluating (Wcr|V^|Wcr)p, after arriving 
at expression (44), we must not perform the sums over a and a', both indices being equal 
to one if spins are assumed to point up. Hence, we are left with quantities ri^ ^ — 
rf\ = 6s',i and the sums over s and s' simply yield 1. Consequently, in Eq. (51), we must 
omit factors 4 and 2 in front of the two sums. Moreover, the cell number instead of N/2 
becomes and the lattice parameter a„, present in Eqs. (52) and (53), becomes ap. For 
these reasons, numerical factors Cd,u and Ce,u must be multiplied by 2 ■ 2^/'^/4 and 2 ■ 2^/'^/2, 
respectively. In conclusion, the energy per particle in the polarized case is 

2 

CpK, r,] =(W„,p|i^|We,,p) / (iV^) (69) 

1 



with 



r, 2 



6\§ l/3\i /6\3 



IV. FUNCTIONAL EQUATION FOR 

From Eqs. (65) and (69) it appears evident that the general form of the energy per particle 

is 

eh,r,] = ^ + ^^ + ^^, (71) 

where we must append to coefficients q, Cd and Ce further index u or p depending on the 
polarization of the system. In principle the function ui is determined by the condition that 
e[i^i] be minimum. For completeness, we derive the equations that follow from this condition. 
First we observe that z/(^) can be written as 



^i(0 = \/ni{Oe'^'^^\ (72) 
where a;(^) is a real function. Then, i:[z/i] can be written as 

t[nuu;]= [ K(f)(Va;)-(Va;) + (Vv^)-(Vv^)]rf=^e, (73) 
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— * — * 

while Vd[vi\ = Vci[n\\ and vJ^U]] = Ve[ni] since these do not depend on uj(^). Functions u;((^) 
and ni{^), with ni(.^) normahzed, are obtained minimizing the functional 

S[n,,uj,X] = At[ui] + Bvdiui] + C Velui] (74) 

where we have put A = Ct/rs"^, B = Cd/rg and C = Ce/r^. Thus, ni(^), uj{^) and A are the 
solutions of the following equations 

0, -^ = 0, (75) 



^ = [ ni(eVe-i = o. 



Observing that 



= -2[n,mM0 + (Vm) . (V^)], (76) 
6t[ni,uj] 



the first two equations of (75) respectively convert into 



(Vc^) • (Vuj) - nr'/^VW^^ (77) 
= r^l(e' + + ni(e'-0, (78) 



and 



n^iOVMO + (Vni(O) ■ (V^(0) = 0, (79) 

f ^ I m 1 2 



>/2yn l£'l 



= Anii/2(a. (80) 

Thus, the SLQCS of form (19) or (68), which yields the minimum expectation value of H, 
is determined by the real functions ^[^^ and n\[^ that solve Eqs. (79) and (80) within 
the cell Vq under the further conditions that u;(.^) and n\[^ are periodic and that ?^i(0 is 
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non-negative and normalized. The periodicity condition can explicitly be accounted for by 
expanding uj{^) and in terms of the complete orthonormal set of functions ^mlO = 

e^"^'^/^/8 and substituting the expansions in the above two equations. We do not report 
the resulting equations for the coefficients of the expansions because the resulting equations 
look more involved than Eqs. (79) and (80). 

The derivation of Eqs. (79) and (80) coincides with that of the HF equations reported in §2 
of [^^] because in both cases one minimizes the expectation value of H over a determinantal 
wave function. However, in contrast with [^^], the SLQCS approach assumes that the 
eigenfunctions are now the functions defined as being equal to z/(^ + m) inside the mth 
cell and zero elsewhere. This fact authomatically ensures their orthogonality and explains 
why Eq.s (79) and (80) involve a single function [z^i(0] ^nd a single eigenvalue [A]. These 
considerations deserve a further remark. One would expect that (Wcr|V^|Wcr) contains a 
direct and an exchange term because |Wcr) is a completely antisymmetric state, but the fact 
that two different eigenfuntions have supports with a void intersection implies that SLQCSs 
by construction yield exchange and correlation contributions equal to zero. 

To better understand the assumptions underlying the SLQCS choice, we observe that a 
normalized quantum state relevant to particles confined to box V can generally be written 
as 

[ ^(^1' ■__ ^^) ^t(^^) . . . ^t(^^)^^^ . . . ^^^|o^^ (81) 

where </'(xi, . . . , x^r) is a normalized and completely antisymmetric function. (For notational 
simplicity, we shall confine ourselves to the full polarized case till the end of this section.) The 
fundamental state is determined by the function </'o(xi, . . . , xjy), implicitly dependent on r^, 
that makes the expectation of H over the resulting |\E'o) minimum. In its fundamental state, 
the system, characterized by a translation invariant Hamiltonian, will have a crystalline 
behaviour at a given value if 00 (xi, . . . , xa?) is such that 

n(r) = (^oln(r)l^o) = (82) 

N I |0o(r,X2, ...,Xiv)prfX2---rfXAr 

7y(8{iV-l) 

turns out to be a periodic function. In this case, state (81) with = 0o is the exact 
quantum crystalline fundamental state and the resulting energy will have a null correlation 
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contribution. It is straightforward to check that states of form (68) correspond to the choice 

(/.o(xi,X2,...,x;v) = 5Z(-l)''x (83) 

p 

(z/(x,J^(xi, G Vi)) ■ ■ ■ (z/(xi^)^(x,^ e Vn)), 

where the sum runs over all possible permutation of {1, . . . , A^} and (— )^ is the parity of 
permutation {ii, . . . ,1^}. According to this relation, function 0o is equal to zero if two 
or more of its variables fall within the same cell. Consequently, by choosing a SLQCS of 
form (68) one implicitly assumes that the pair correlation function cr,p\n{vi)n{T2)\^ cr,p) i 
[i.e. the probability density of finding two particles at two different points ri and r2) is 
respectively equal to zero or n{yi)n{v2) depending on whether ri and V2 lie in the same or in 
different cells. This behaviour is radically different from that obtained by (81) that neither 
factorizes nor, more importantly, vanishes if ri and Y2 lie in the same cell. For this reason, 
the SLQCS approach fails in describing the "fluctuations" in the electron density. One 
expects that when these are important the SLQCS results will not be accurate. However, in 
the limit of a very narrow n(r), the SLQCS aforesaid beahaviour is similar to that of an ideal 
classical crystal with a point-like particle within its unit cell so that the SLQCS description 
can be realistic only within the region where the system is crystalline and nir) turns out 
to be rather narrow. 

V. A NUMERICAL APPLICATION 

Instead of numerically solving Eqs. (79) and (80), which appears to be a rather involved 
problem, we shall determine the region where inequality 

(^,,|H|^,,) < (F„|H|F„) (84) 

is obeyed confining ourselves to SLQCS defined by a particular class of functions. 
Clearly, the resulting region will be smaller and the expectation value larger than those 
obtained by solving Eqs. (79) and (80). Nonetheless, the results appear to be interesting. 
Hence, we shall assume that z^i(C) has the functional form 

z/i(e, a) = G(ei, a)G(6, a)G{^s, a) (85) 
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with 



C{a) 



v/^^(-,0;2«^ 



-1/2 



(86) 
(87) 



where a is a real positive parameter to be determined by a minimization procedure and 
0; la) is a speciahzation of \I/(a, c; 2;), the confluent Hypergeometric function deflned in 
§2.5 of Ref. p^]. The reported C[a) expression ensures that 



(88) 



Functional choice (85) implicitly assumes that = 0. Further, we assume that i'i{C,,a) 
is an even non-negative function, fully symmetric in variables ^1,^25^3, the components of 
C, varying within the interval [—1, 1]. These assumptions imply that we are considering a 
crystal with a simple cubic structure. Moreover, functional choice (86) implies that z/i(^,a) 
and all its partial derivatives approach zero as ^ approaches the cell boundary. Then, their 
FTs decrease rather quickly as the momentum increases and the series involving the FTs 
can safely be truncated. The factorized dependence of ui on ^1,^2 and ^3 simplifles some of 
our results. In particular it results that 



t(a) 



6 



/[ 

^0 



91 



(f,«)=n^i(o, 



a 



6 



Trrrij, a) 



,n5=i[G'2(vrm„a)]^ 



m 



(89) 
(90) 

(91) 

(92) 



v,{a) = -48 / Gii^i,a)d^i 
Jo 



€1 



Gi(^2,a)c?6x 



veF+er+ei 



where 



G2(g,a) = 2 / cos{qOG\i,a)di. 



(93) 

(94) 
(95) 
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FIG. 3: Left: unpolarized case. Continuous curve: Log\Q\ct^yt{a)\ vs. LogiQ^a); short dash: 
Cd,uVdi<^)j dot-dash: Ce,«t'e(a); long dash: Cd^uVdio:) + Ce,uVe{o). Right: the same quantities for the 
polarized case; the horizontal axis values must be shifted by -3. Finally, the two dotted curves 
show the best a value [as LogiQ{a) on the horizontal axis, in terms of LogiQ{rs), on the vertical 
axis] for the unpolarized (left) and polarized (right) case, respectively. 

The above equations makes it evident that the determination of t{a), Ve{a) and Vd{<y) is nu- 
merically straightforward. In fact, the integrals must be performed over compact sets where 
the integrands are continuous while series (91) can safely be truncated because G^(g, a), as 
shown in Appendix A, asymptotically decreases as 

vl/(i/2,0,2a)g3/4 '^^^^ " + ^r/S). (96) 

In practice, integrals (89), (94) and (95) have been numerically determined considering 10^ 
intermediate points; the ^i, ..,^3 spacing in Eq. (93) was taken equal to 10^'' and the largest 
absolute value of each component of m in Eq. (92) equal to 30. Moreover, quantities t{a), 
Vd{<y) and Ve{a) were evaluated over a grid of a values that span the intervals (0, 1] and 
[1, 10] with spacings respectively equal to 0.025 and 0.25. The correspondent values are 
denoted by taj = ta{aj), Vd,a,j and t'e.aj, with a = u,p and j = 1, . . . , 80, and are plotted in 
Fig. 3. Once they are substituted in Eqs. (65) and (70) one respectively obtains quantities 
eu{aj,rs) and ep{aj,rs)- Finally, the energies per particle in the unpolarized and polarized 
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case respectively are 

e«;cr,n(rs) = minje„(aj,r^) (97) 

and 

ewcrA^s) = mmjep{aj, r^). (98) 

The two minimizations also determine the a value associated to a given value in the 
unpolarized and polarized case. The corresponding a„(rs) and a{rs) functions are plotted 
in Fig. 3. The quantities eu,cr,«(rs) and ewcr,p{rs) are plotted in Figs. 1 and 2 together with the 
other approximations mentioned in sect. I. Their comparison with the relevant HF energies 
shows that the SLQCS approach predicts that the unpolarized and polarized jellium models 
respectively acquire a simple cubic structure as becomes greater than 90 and 28, and that 
the stable phase is the polarized one as Vg > 8.5. The most recent QMC result^^ locates 
the jellium model's transition from the homogeneous phase to the polarized bcc structure 
at Vs = 106 ± 1. Further at = 100, the energy per particle [in Rydberg] is: -0.0112 by 
Eq. (3); -0.01448 by the SLQCS approach; -0.014986 by the general HF approach^o, -0.01526 
by Eq. (9) and -0.0153530(8) by QMC^^. These values show that the relevant correlation 
contributions decrease, i.e. electron interactions are better accounted for, as we go from the 
first to the last approximation. Moreover, the fact that the values of the second and third 
approximation are nearly identical indicates that the SLQCS is close to the crystalline HF 
one, which also predicts a vanishing n{r) on the cell border (see Fig. 5 of Ref.[^^]). A further 
comparison is possible with the bcc QMC values^^ at = 125, 150 and 200 (in the last case, 
we consider the energy values resulting from Eq. (12) of Drummond et al.'s). In the three 
cases the relative errors of the SLQCS values with respect to the QMC ones respectively 
are: 4.9%, 4.3% and 2.0%. These results show that as the dilution increases the SLQCS 
more completely accounts for all electron interactions. The reason is probably the following. 
The dotted ap(r<,) curve, reported in Fig. 3, shows that the sharpness of the resulting particle 
density increases as Vg increases and this fact, as explained at the end of the previous section, 
makes the SLQCS approach more reliable. 

VI. CONCLUSIONS 

Similarly to the HF analysis of^*^, this paper provides another application of the Rayleigh- 
Ritz principle to the Wigner crystal by considering the exact many-body Hamiltonian, in 
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contrast with Wigner's and Carr's analyses. To this aim, we considered a class of quantum 
crystalline states that are essentially determined by the electron density within the unit cell 
and, by construction, allow no fluctuation in the electron number within each cell. These 
SLQCSs have the advantage of making the numerical evaluation of the expectation value 
of the Hamiltonian very simple if one assumes a convenient functional form for the electron 
density. In this way, one easily obtains an upper bound for the energy of the fundamental 
state of the jellium model. The comparison of this bound with that obtained by the HF 
procedure, also based on the variational Rayleigh-Ritz principle, shows that the HF bound is 
better for the reasons explained in sect. IV. However, the comparison of the SLQCS energies 
with the HF and the QMC ones shows that the SLQCS values approach to the latter as the 
density decreases and this indicates that the fluctuations in the electron density within each 
cell decrease with the density. Finally, the fact the SLQCS approach is a HF one implies 
that its results can be improved by the perturbative procedure of MoUer and Posselt^^. 
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APPENDIX A: ASYMPTOTIC BEHAVIOUR OF CP{q,a) 

To show that, as g — cx), Eq. (96) is the leading asymptotic term of G'^{q, a), defined by 
Eq. (95), we put 



Denoting the above integrand by f{t, a), one sees that / is an analytic function throughout 
the complex plane t = u + iv deprived of the points 0, 2 and oo. Integral (A2) can be 
written as 



g{q,a) = G^{q,a)/C^{a) = 2e 



3fJ[e-^'?F(g, 2a)] 



(Al) 



with 




(A2) 




l+i 



(A3) 



+ t 




/(I + iv, a)dv. 
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In the limit a oo, one easily verifies that the second integral vanishes because the integral 
exponentially decreases if g > 0. The integrand of the third integral reads 

SO that, in the limit a — > oo one finds 



This expression, multiplied by e as required by Eq. (Al), is an imaginary quantity and, 
consequently, it does not contribute to g{q,a). Thus, we find that 

POO 

g{q,a) = 2e^'^^[e-''^ / h{v,a)dv\, (A4) 

As g ^ oo, the only region v ^ contributes to integral (A4), so that the integrand can fairly 
be approximated by e-"/4 + «a/2i' g^j^j ^j^g ^-j^g asymptotic behaviour of g{q,a) is determined 
by that of 

/•oo 

J{q,a)= / e-^^+'^/^-^iy (A6) 
Jo 

due to the relation 

g{q, a) ^ 26^"'^ [e-''e-"/^J{q, a)] . (AT) 
Considering the new integration variable t = v^fqfla we find that 

xJxix) _^fqJ{q,(^) 



oo 

e~ 







2a .J2a 

e-"'^(*))cit (A8) 





where we have put = {t — i/t) and x = y/2aq. We apply now the saddle point method 
(see, e.g., Ref. f^]). Thus, we look for the points of the complex plane t = u + iv where the 
derivative of h{t) vanishes. These points are 

=e— /4 = i^ and = e^'""/' = (A9) 
V 2 V 2 

The steepest descents through the above two points are determined by the condition 
= const. In particular, the two steepest descents through ti respectively have para- 
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metric equations 



uAv) = 0<v < -I/V2, (AlO) 



= i±^, -l/V2<v<-V2, (All) 
2{v + v2) 

with 

A{v) = 1 - 8v'^ - 8V2v^ - iv^. (A12) 

Moreover they form a continuous contour, the first descent going from the origin (when 
V = 0) to ti (when v = —I/V2) and the second from ti to {00 — i\pi) (when v = —a/2). 
Considering the further hnear contour joining (00 — to (00 + i 0), integral (A8) can be 
written as 



e 







+ / e-"'^("^(")+'")((u;(tO + z))rfi;+ 

J-V2 

The last integral vanishes because x > and h{oo + iv) = 00 in the considered integration 
domain. Further, as w — and as f ^ — "\/2 it results that '^[h{ui{v) + iv)] — > cxo and 
3?[/i(M2(f ) + iv)] 00. Since we are interested in the large x behaviour, the main contri- 
butions to the remaining two integrals in Eq. (A13) arise from the f-values very close to 
— 1/\/2. Expanding 'Si{h{ui{v) + iv)) and 'St{h{u2{v) + iv)) around v = — 1/-\/2 one finds 

^h{ui{v) + iv)) =^{h{u2{v) + iv)) 

V2 + 2(2 + ^2) {v + l/V2f + o 

and integral (A8), equal to Eq. (A13), asymptotically behaves as 



/oo 
-00 



(1 + 72-,.— Ov-^^-^e-^'. ,A14, 

Combining Eqs. (A14), (A8), (A7) and (Al), one finds that the leading asymptotic term of 
G'^{q,a) is that specified by Eq. (96). 
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